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Abstract 

The leading-log QCD corrections to K — K mixing in R-symmetric supersymmetric models 
are computed using effective field theory techniques. The spectrum topology where the gluino is 
significantly heavier than the squarks is motivated and focused on. It is found that, like in the 
MSSM, QCD corrections can tighten the kaon mass difference bound by roughly a factor of three. 
CP violation is also briefly considered, where QCD corrections can constrain phases to be as much 
as a factor of ten smaller than the uncorrected value. 
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1. INTRODUCTION 



It has recently been proposed in pQ that by extending the usual minimal supersymmetric 
standard model (MSSM) to a theory with an additional U(1)r symmetry, the SUSY flavor 
problem can be solved, at least in part. In particular, the additional R symmetry forbids 
several flavor-changing neutral current (FCNC) operators usually found in the low energy 
effective theory of the MSSM; and those operators that are still there can be suppressed in 
large regions of parameter space, even for sizable flavor mixing angles. This is in contrast 
to the usual requirement of a "flavor-blind" mechanism, such as gauge mediation, to forbid 
the presence of such mixing angles. 

Models with a (partial) U(1)r symmetry were considered several years back in [2j (3J H], 
where the gauginos were given a Dirac mass. Several papers have extended on this idea 
over the years |5l El El El El HDl [HI H21 H31 HH H51 US]- What was shown in was that in 
regions of parameter space where the gluino is much heavier than the squarks, the FCNCs 
are under control even for sizable mixing angles. Fortunately, it seems that such spectrum 
topologies are rather generic in models with Dirac gauginos [5]. 

Of course, a viable model of supersymmetry cannot have a completely unbroken R sym- 
metry, since such a symmetry is unavoidably broken by the gravitino mass which arises due 
to the Super-Higgs mechanism whenever you break supersymmetry. This will inevitably 
generate R symmetry breaking in the low-energy sector via the Giudice-Masiero mechanism 
or anomaly mediation or some such device. However, one can imagine that such effects 
are somehow hidden or rendered irrelevant by an as-yet-unknown UV completion. Some 
thoughts along these lines were presented in [TJ. 

Even with these criticisms, the prospect of being able to solve the SUSY flavor problem 
without resorting to the usual flavor-blind requirements is exciting, to say the least. In pQ, 
a tree-level analysis of flavor changing effects from gluino box diagrams was considered. In 
this note, we compute the one-loop QCD corrections to these results. This is believed to 
be the most important correction to K — K mixing, although there are scenarios where 
electroweak box diagrams can also be important. We will follow closely the outline of [17J, 
where a similar analysis was done for the usual MSSM. 

2. THE EFFECTIVE LAGRANGIAN 

The R symmetry can be seen to give some very powerful constraints on the low-energy 
MSSM: 

• No Majorana masses for the gauginos. 

• No A-terms among the scalar fields. 

• No \x term for the Higgs, although there is a term. 

Putting these facts together leads right away to the result that the only coupling between 
the "left-handed" squarks and the "right-handed" squarks (and sleptons) is through the 
usual Yukawa interactions, and these are flavor diagonal. Therefore, any flavor-changing 
couplings can only come from the soft squark mass matrices, and these are necessarily of 
the form "left-left" or "right-right". In what follows, we will assume that the squarks are 
nearly degenerate in mass (this is a reasonable assumption for the first two generations in 
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an mSUGRA-like scenario), and parametrize the off-diagonal elements of the mass matrix 
in terms of dimensionless numbers 1 : 



where Q is the left-handed squark and d is the right-handed down-type squark and Mq is the 
universal squark mass. These are the only relevant squark-mixing angles for K — K mixing. 

To compute the amplitudes in a process with several scales, the best way to proceed 
is through the construction of an effective theory, where we construct operators of lowest 
possible dimension using only the low energy degrees of freedom, in this case the quarks. To 
this end, we can write down an effective Lagrangian of the form 

^ 9 ' n 

where the leading operators 2 are dimension 6: 

Oi = (47 M 4) (4 7 ,4) (3) 

o A = (44) (44) (4) 
a = (44) (44) (5) 

i, j are color indices, and there is a similar operator 0\ with L <-> R. We are using the same 
numbering scheme and normalization as [T7j which was also used in [I]. Notice that the lack 
of any "left-right" couplings means that 62,3 = in their operator list. The /i in Equation [2] 
is the renormalization scale, the dependence of which cancels between the operator and the 
Wilson coefficient. This implies that the Wilson coefficients must obey the Renormalization 
Group Equation: 

dC m rp 

^~dy = 7mn o 

where 7 mn = (■^~ 1 )ma ^~" is the anomalous dimension matrix. 

There are three spectrum topologies that one can imagine: (1) Mq ~ Mg, (2) Mq ^> Mg 
and (3) Mq <C Mg. The first two cases turn out to be identical to the usual MSSM done 
in [T71 US], with 5lr = f§(x) = 0. However, case (3) turns out to be quite different. Since 
these models are most favorable in this case, we will present this calculation in detail. 

We proceed in four stages. The first step is to integrate out the gluinos at the scale 
M^, where we match to an intermediate theory with quarks and squarks as the degrees of 
freedom. We will then run the theory down to the scale Mq, which allows us to resum 
the leading-log QCD corrections between the two scales. Next, we match our intermediate 
theory to Equation [2} Finally, we can do the usual QCD running of this theory down to the 
hadronic scale, integrating out the heavy quarks as we go. Ideally we would like to run the 
theory down to the mass of the K meson, but QCD is strongly coupled at such scales, so we 
will stop the running at a scale /i had ~ 1 GeV defined by the condition that a s (/Xhad) = 1- 



1 These quantities are usually referred to as SlLtSrr in the literature, but since there is no flavor mixing 
of the "left-right" form, we drop the second L(R) for notational simplicity. 

2 In this list the fermions are 4-component fields, whereas in the rest of this paper, we use 2-componcnt 
notation. It is very easy to go from one to the other: just replace ipR i>R and leave the ipL fields alone. 
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FIG. 1: Integrating out the gluino. In this and all Feynman diagrams, momentum is flowing left 
to right, and the arrows on the fermion lines represent spinor flow. 



2.1. Step 1: Matching to the Intermediate Theory 

2.1.1. General Operator Analysis 

The gluino couplings to the quark-squark are given in the full theory by 

AC = -V2g s [(q*T a q)\ a + \ a (qT a q)} (7) 

Here T a are the 577(3) generators in the fundamental representation, and A a is the Weyl 
spinor field representing the part of the gluino that interacts directly with matter. Since it 
no longer has a Majorana mass, the only propagators are of the form (AA) and the diagram 
in Figure [T] gives 

-id ■ (pi 



iM 



c% 2rparpa 
^xs ij mn 



P3) 

-Ml 



(q?*ql 



(8) 



This amplitude then matches onto an intermediate effective theory with six operators up to 
flavor structure: 

9 2 s(M- 9 



Ml 



(9) 



where 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

Here, P M is the QCD covariant derivative required by gauge invariance, and each operator 
has a Wilson coefficient: 



Qi = 








0,2 = 


[(iV^qi] 




[qi^4]((W^yqf) 


Q 3 = 






iq\^(W,q 2 y}(qiq?) 


Q 4 = 




W 2 ](f 3 C) + 


[qt^(iV„q 2 y}(<f 3 qt) 


Q 5 = 






[^qiMiW.qiy) 


Q 6 = 






[^qiMiiV.qiy) 



D\{Mg) = -l + 0(g 2 s ) (16) 

Ai(M s ) = +l + 0(g 2 s ) (17) 

Ds-i{M s ) = 0(g 2 s ) (18) 

Note that unlike the situation in [T7j these operators are dimension 6, not dimension 
5. This point was emphasized in [T| and stems from the key fact that there is no gluino 
Majorana mass, and hence no (AA) propogator. 
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2.1.2. Class A versus Class B Operators 



In order to match the intermediate theory to Equation [2] we must specify the flavor 
structure of the operators. There are four flavor combinations that we will need, and these 
can be divided into two classes of operators depending on their color structure. The "Class 
A" operators corresponding to Equation [TO] are: 



q\ = [(^)^4](4€) + [*L^£){(a>vd L ys£) (19) 

Ql = [(iP^)W K ](44*) + [4^4]((zP M 4)^4*) (20) 

and the "Class B" operators are: 

Q? = + [d\a^]((iV,d L yd%) (21) 

Q\ = [(^ R y^si](~s j R st) + [s R ^si]((iD,~s R y~st) (22) 

and similarly for Qi^Q. Since QCD is vector- like and flavor blind, all of these operators will 
share the same matching coefficients Di(Mg). However, since the color flow is different 3 , 
the Class A operators will run differently than the Class B operators: see the appendix for 
details. 



2.2. Step 2: Running the Intermediate Theory 



To compute the one-loop running of the operators we must evaluate the diagrams 
shown in Figure [6j Performing this calculation is a straightforward (if tedious) exercise and 
yields the result for the anomalous dimension matrix presented in the appendix. 

Now we can compute the Wilson coefficients at the squark mass scale, using bo = — 5 



n 



-77/5 



where we define 



V 



(23) 



(24) 



OC S {Mg) 

and 7 is the anomalous dimension matrix without the a s /2n prefactor. The index I is for 
whether the operator is a Class A or a Class B operator. The solutions are quite messy and 
are computed numerically in what follows. 



2.3. Step 3: Matching to the Low Energy Effective Theory 

We can now use the renormalization group to evolve the operators from the gluino mass 
scale down to the squark mass scale, and then consider the diagram where two operator 
insertions and two flavor-violating couplings generate the dimension 6 operators in Equation 



3 This is only true for 21,3,5; for 22,4,6 the color flow is identical between Class A and Class B operators. 
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FIG. 3: Matching calculation to involving Class B operators. 
[2j see Figure |2](3j This matching will proceed as: 

01 a 

C 1 (M,)(0 1 (M,)) = ^ E (D?(Me)D}{M 9 )) (Q\(Mq)Q)(Mq)) (25) 

21 fi 

(7 1 (M,-)(a i (M 9 -)) = — 2 (Df(M 9 )D}{M 9 )) (Q 2 (M q )Q 2 (Mq)) (26) 

C 4 , 5 (M ? -)(0 4j5 (M 9 -)) = ^ E iPfWBf ( M .-)) (2, 3 (^-)Q, 4 (M 9 -)) (27) 

Since the do not contain derivatives we can safely match matrix elements involving states 
of quarks with vanishing momenta, where we immediately see that (^3,4) = 0. So these 
operators do not contribute in the matching, although they do contribute to the running. 
By doing the loop we calculate the low energy Wilson coefficients at fi = Mq to be 

d(M 3 -) = -9S 2 L ((3D 2 + 2D 1 D 2 + Dl + 3D 2 5 +2D 5 D 6 + D 2 ) 

-2{3D 1 D 5 + DiD 6 + D 2 D 5 + D 2 D 6 )) \^ =M _ (28) 

C^Mq) = -72 S* L 8 R {D 1 D 2 + D 5 D 6 - D 1 D 6 - D 2 D 5 ) | M=M _ (29) 

C 5 (Mq) = -36 S* L S R ((D 2 + D\ + D\ + D\) - 2(D 1 D 5 + D 2 D 6 ))\^ =M _ (30) 

and C\ = C\ with Sl <-> Sr. Here it is understood that the Wilson coefficients in Equation 



28 are Class A while the Wilson coefficients in Equations 29 30 are Class B. Plugging in the 



uncorrected values for the Wilson coefficients in Equation 16 - 18 gives agreement with [TJ. 
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m s ttik 


fK 


Arrix 


Bi 


B 2 


#3 


7 MeV 


95 MeV 498 MeV 


160 MeV 


3.48 x 10" 12 MeV 


0.6 


1.03 


0.73 



TABLE I: Numerical values for the parameters in the matrix elements [19] . 



2.4. Step 4: Running to the Hadronic Scale 



Once we have the Wilson coefficients at the squark scale, we can run the scale down to 
A*had- To do this, we must calculate the anomalous dimensions for these operators. This has 
already been done in [IT] : 



7i 
"Hi 



71 

«s / -8 
2vr V -3 +1 



i,j = 4,5 



(31) 
(32) 



Along the way we integrate out the quarks. This will not generate any new operators at 
this order, but it will change the running of a s . 

Finally we must calcuate the remaining matrix elements at the hadronic scale. They are 
given by 



(K\0,(fx had )\K) 



(K\0 5 (fx had )\K) 



1 1 

2 + 3 



m K 



m s + m d 
m K 



mKfK B A 



m s + m d/ 

where Bi are the bag factors characterizing the vacuum saturation approximation. 



(33) 
(34) 

(35) 



3. RESULTS 

3.1. K — K Mass Difference 

Now that we have computed the leading contributions to the effective Lagrangian in 
Equation [2] we can proceed to place bounds on the model by looking at the K — K mass 
difference given by 

Am K = 2Re ((K\C eS \K)) (36) 

This number is experimentally constrained to be less than 3.48 x 10 -12 MeV. In what 
follows we will use the numerical values for the matrix element parameters given in Table |T} 
In addition, we use a s {Mz) = 0.118. 

In Figure [4]{5] we plot the upper bounds in the Mg — Mg plane for various values of 5l 
and 5r. In Figure [4] we show the plot for a universal value of 5. It is very clear from these 
plots that the QCD corrections are very significant, roughly a factor of three stronger. This 
is similar to the results of [17J. This somewhat weakens the resolution to the flavor puzzle 
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FIG. 4: Exclusion plots in the Mg — Mq plane. The dashed lines are the uncorrected bounds, while 
the solid lines include the QCD corrections. The region below the lines is excluded. Going left to 
right: 5 L = 5 R = 0.03, 0.1, 0.3, 1.0. 



in [T], but it does not destroy it. In particular, for 5 ~ 1 the model is strongly constrained, 
but even for 5 ~ 0.1 it is still very open. 

Of course, for 5 ~ 0(1) our approximations formally break down, since such a large off- 
diagonal mass term would imply a larger splitting between the d and s squarks. While we 
do not expect this to change the bounds significantly, it is worth remembering when trying 
to do a careful phenomenological analysis of the model. 

In Figure [5] we do the same analysis but with 5r = 0. In this case only C\ is nonzero. As 
is clear from the figure, the QCD corrections do very little to the constraints. 

That we roughly agree with [IT] is not unexpected. In both cases, it is the low energy 
physics below the squark scale that dominates in the running. You can guess that this 
is the case since at lower scales (1) a s is larger and (2) the resummed logs are larger: 
log(m^/mg) log(mg//^had) [20]. Nevertheless, is is good to show the full calculation to 
give quantitative results. 
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FIG. 5: Same as Figure[|J but with S R = 0. From left to right: 5 L = 0.1,0.3, 1.0. 



3.2. CP Violation 



In addition, there is also a nontrivial phase that can enter through the product 5^5 r, so 
we can consider the effect of this phase on the CP-violating parameter 

|£Kl = V2A,n K (37) 

We can read off the bounds on CP violating phases from Figure [1] as explained in [Tj . If 
we let 6 = arg(52^j?) <C 1, we have that ek oc 6. If we use the example points in [1], 
M- g = 3500 GeV, M~ q = 400 GeV, and 5 L = 5 R = 0.06 (0.25), and assume that the SUSY 
contribution saturates the bound 4 , we find that 9 < 9.8 (0.57) x 10~ 3 . This is an order of 
magnitude stronger than the uncorrected results and suggests that the phase might have to 
be unnaturally small in order to avoid experimental constraints. 



4. CONCLUSION 

In this paper, we have computed the leading-log QCD corrections to the model discussed 
in pQ. It was realized then that an R— symmetry built into a SUSY model with Dirac 
gauginos naturally avoids the flavor puzzle found in typical SUSY models. We find that 
QCD corrections can significantly affect this result by nearly a factor of 3, similar to the 
usual SUSY story in [17]. The model is still significantly less fine tuned, with 8 ~ 0.1 
allowed in a typical models, although in the universal case, CP violation can still be very 
constraining. It is clear that QCD corrections are important in general and must be included 
when studying FCNC processes in R-symmetric supersymmetric models. 



APPENDIX A: ANOMALOUS DIMENSION MATRIX IN THE INTERMEDI- 
ATE THEORY 

The intermediate theory is described up to flavor structure by the six dimension-6 op- 



erators given in Equations 10 15 To compute the QCD anomalous dimension matrix, one 



This is what was done in [T|. 
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FIG. 6: QCD loops that contribute to the leading order counterterms to Class A operators. The 
Class B operators are similar. 

must evaluate the diagrams in Figure [6] for each of these operators. The Class A operators 
<2^3 5 have a different color structure than the corresponding Class B operators Qi' 35 and 
will therefore run differently. However, the operators Q^^e' have the same color structure 
and therefore both classes run the same way; this has been confirmed by direct calculation. 
The final result is a 6 x 6 matrix for each class: 
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39 
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3 
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+ 15 
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-n 
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14 
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-6 


+2 


10 
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13 
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91 
12 
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" t_ 12 
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91 
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(Al) 
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